import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D

vertex = np.array([[0, 0, -0.371091],
                   [0, -5, -0.323287],
                   [0, -10, -0.275484],
                   [10, 0, -0.252012],
                   [10, -5, -0.243377],
                   [10, -10, -0.234741],
                   [20, 0, -0.0725019],
                   [20, -5, -0.0945335],
                   [20, -10, -0.116565],
                   [30, 0, -0.01207],
                   [30, -5, -0.025601],
                   [30, -10, -0.0391319]])

gauss = np.array([[7.88675, -3.94338, -0.151576],
                  [7.88675, -1.05662, -0.251061],
                  [2.11325, -1.05662, -0.429909],
                  [2.11325, -3.94338, -0.330425],
                  [7.88675, -8.94338, -0.237886],
                  [7.88675, -6.05662, -0.142655],
                  [2.11325, -6.05662, -0.285558],
                  [2.11325, -8.94338, -0.380789],
                  [17.8867, -3.94338, -0.0643234],
                  [17.8868, -1.05662, 0.000276597],
                  [12.1132, -1.05662, -0.151318],
                  [12.1132, -3.94338, -0.215918],
                  [17.8867, -8.94338, -0.0480461],
                  [17.8868, -6.05662, -0.0944682],
                  [12.1132, -6.05662, -0.233106],
                  [12.1132, -8.94338, -0.186684],
                  [27.8868, -3.94338, -0.0610773],
                  [27.8868, -1.05662, -0.026193],
                  [22.1132, -1.05662, 0.00106147],
                  [22.1132, -3.94338, -0.0338228],
                  [27.8868, -8.94338, -0.00400123],
                  [27.8868, -6.05662, -0.0528104],
                  [22.1132, -6.05662, -0.0485451],
                  [22.1132, -8.94338, 0.00026404]])


fig = plt.figure('3D scatter plot')
ax = fig.add_subplot(111, projection='3d')

ax.scatter(vertex[:, 0], vertex[:, 1], vertex[:, 2], c='r', marker='o')
ax.scatter(gauss[:, 0], gauss[:, 1], gauss[:, 2], c='g', marker='^')

plt.show()

